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ABSTRACT 



Context. We present the CoRoT light curve of the/? Cep star HD 180642, assembled during the first long run of the space mission, as 
well as archival single-band photometry. 

Aims. Our goal is to analyse the detailed behaviour present in the light curve and interpret it in terms of excited mode frequencies. 
Methods. After describing the noise properties in detail, we use various time series analysis and fitting techniques to model the CoRoT 
light curve, for various physical assumptions. We apply statistical goodness-of-fit criteria that allow us to select the most appropriate 
physical model fit to the data. 

Results. We conclude that the light curve model based on nonlinear resonant frequency and phase locking provides the best represen- 
tation of the data. The interpretation of the residuals is dependent on the chosen physical model used to prewhiten the data. 
Conclusions. Our observational results constitute a fruitful starting point for detailed seismic stellar modelling of this large-amplitude 
and evolved /3 Cep star. 
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1. Introduction 

The B1.5II-III star HD 180642 (variable star name V1449 Aql, 
Vmag - 8.29) was ide ntified as a candidate new /3 Cep star by 
IWaelkens et alj (1 19981) f r om H ipparcos data. This classification 
was confirmed by lAertsI (12000'). who identified the detected fre- 
quency of 5.4871 d~^ (63.508/vHz) as a radial mode with a large 
amplitude of 39 mmag in the V band, by interpreting amplitude 
ratios derived from multicolour Geneva photometry obtained 
with the P7 photomultiplier instrument attached to the 0.70m 
Swiss telescope at La Silla, Chile. 

Given that HD 180642 is the only known /3 Cep star with 
appropriate magnitude in the field -of-view of the CoRo T space 
mission (see "The CoRoT Book": iFridlund et alj|2006h . we un- 
dertook a preparatory observing effort to assemble data to be 
added to the CoRoT light curve. Several high-resolution spectra 
were taken with the FEROS@2.2-m ESO/MPl telescope in La 
Silla, Chile, in 2005. The latter led to an estimate of the funda- 
mental parameters of the star: T^ff = 24 500 ± 1 OOOK, log g - 
3.45 + 0.15, as well as an overall line broadening of 44kms"' 



* The CoRoT space mission was developed and is operated by the 
French space agency CNES, with participation of ESA's RSSD and 
Science Programmes, Austria, Belgium, Brazil, Germany, and Spain. 
Send offprint requests to: P. Degroote 

** Postdoctoral Fellow of the Fund for Scientific Research, Flanders 
Send offprint requests to: P. Degroote 



and some evidence for a mild nitrogen excess dMorel & AertsI 
|2007), as disc overed in o ther ySCep stars from high-precision 
spectroscopy (iMorel et al.l l2006, 2008), which could betray the 
existence of deep mixing. The combination of the low gravity 
and high pulsatio nal amplitude of th i s clas s member is rather 
exceptional (see iStankov & Handled ( l2005b . their Fig. 8) and 
seems to suggest an object near the end of the core-hydrogen 
burning phase, almost ready to cross the Hertzsprung gap in the 
Hertzsprung-Russell diagram. 

Adding the space photometry from Hipparcos to CoRoT and 
ground-based data, brings the total timespan of observations to 
18 years. The dominant mode of the star is present in all of these 
datasets and we thus have the means to determine its frequency 
stability over time. On the other hand, the high time-sampling of 
the CoRoT light curve combined with its low noise level, give 
us the possibility to look for variability far beyond this domi- 
nant mode. The richness of the CoRoT frequency spectrum led 
at once to the conclusion that the monope riodicity of the star 
must be refuted, as was akeady suggested bv lUvtterhoeven et al.l 
(.2008i) from the ground-based data. 

We have entered a new stage of precision in observational 
astronomy with the CoRoT space mission. In this paper we thor- 
oughly investigate the variability of HD 180642 from single- 
channel photometry. Additional time series of multicolour pho- 
tometry and high-re solution spectroscop y of the star are the sub- 
ject of a twin paper (iBriquet et al.ll2009 ' 
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Fig. 1. Part of the reduced CoRoT light curve (upper panel), containing 379 785 datapoints in total. Despite the intrinsic equidistant 
nature of space based ob servations, g aps are unavoidably present, mainly due to the regular passage of the satellite through the 
South Atlantic Anomaly dAuvergne et al 2009.) . The lower panel shows the CoRoT light curve of the constant star HD 181072. 



2. Observations 

2.1. The CoRoTdata 

The raw light curve from the CoRoT database contains 422 949 
datapoints, with an average time sampling of 32 s during 
156.6 days and starting on fo = HJD 2454232.091674. This 
brings the Nyquist frequency up to 1350d"'. To obtain the 
highest possible precision, roughly 10% of the datapoints were 
deleted because of flagged datapoints (9.8%) and extreme out- 
liers (0.5% have an estimated error value above the 6cr level), 
keeping 379785 datapoints (Fig. [Til. 

After rigorous tests, we decided not to interpolate the re- 
maining points, because the improvements of the spectral win- 
dow do not weigh up against the introduced uncertainties con- 
nected with the gap filling model. The highest amplitude in the 
window function is only ~8% of the main amplitude (Fig.|2]i. 
This means that we effectively spread out the power of each peak 
over several peaks, mainly well separated by ~ 13 d"'. 

A raw estimate of the noise level of the light curve, computed 
as the average of the periodogram between 30 d"' and 40 d"' is 
at 57 yumag or 0.00536 % in relative flux units, and slowly decays 
at higher frequencies to 24yumag or 0.0026% between 100 d"' 
and llOd We do not convert the light curve to magnitudes 
because the transformation from flux is not uniform and the the- 
oretically predicted variations in first order are only linear in flux 
units. An exception is made in the case where the CoRoT light 
curve is used in combination with ground based observations. 
This does not pose a problem because of the dominant mode's 
large amplitude. 

The final reduced version of the light curve has also been 
corrected for long term trends: among an exponential, parabolic 
and linear trend, the linear trend resulted in the best fit, reduc- 
ing significantly the power in the periodogram at low frequen- 
cies. An instrumental cause of the trend seems most probable, 
although long term (periodical) variations in the brightness of 
the star cannot be excluded from this time series alone. Previous 
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Fig. 2. Spectral window of the CoRoT measurements of 
HD 180642 after removing flagged datapoints and outliers. The 
highest peak is 8% of the main peak, and is well separated from 
it. 



ground based observations disfavour the second possibility, but 
do not exclude it, however. 



2.2. The noise properties of the CoRoTdata 

The CoRoT data products contain information on the standard 
deviation of the star's intensity per second, which is interpreted 
as the noise on the data. For the use and interpretation of the ap- 
plied data analysis tools, it is vital to have a good understanding 
of the noise properties. We divided the error by the local aver- 
age flux value of the light curve represented by a polynomial 
fit. A comparison between the Fourier periodogram of the inten- 
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sity measurements and the derived standard deviations, shows 
that in the case of HD 180642, the latter are contaminated by 
the former, and are thus not a reliable estimate of the instrumen- 
tal noise. Instead, we use the measurements of the constant star 
HD 181072 of spectral type A2 and visual magnitude of 9.14, 
which has been observed simultaneously on the same CCD, as 
an appropriate representation of the noise properties. 

Traditionally, uncorrected homoscedastic white Gaussian 
noise is assumed in frequency analyses of ySCep stars. If the 
number of data points is large enough, this implies that the 
distribution of the normalised Fourier periodogram can be 
reasonably well appro ximated by an exponential distribution 
dSchwarzenberg-Czern v 1998). In the following paragraphs, we 
show that none of the assumptions are strictly true, but the devi- 
ations are so small that the classical methods can still be applied 
provided that a correction for correlated data is used. 

By binning the noise measurements in samples of 1000 
points, we can see that the noise is not uncorrected nor ho- 
moscedastic: we identify a continuously rising trend of (1.71 ± 
0.3) X 10"^ percent d"' and a small temporary bump around day 
130 (Fig. [3]). From a log-log plot of the Scargle periodogram, it 
is apparent that the noise is not white: at very low frequencies 
(< 0.1 d"') there is some power excess due to the correlation ef- 
fect. However, white noise is a good enough approximation for 
/>0.1d-i. 

Drawing random samples of 1000 points, reveals that the 
noise is also not Gaussian: the sample mean is consistently 
higher than the sample median, suggesting that the noise has 
nonzero skewness. When falsifying samples of 1000 consecu- 
tive points for normality, by testing simultaneously for skew- 
ness and excess kurtosis, 65% of the samples are rejected at a 
p = 0.01 acceptance level. Bootstrapping the same number of 
samples of 1000 points yield a rejection rate of 85%. A Gaussian 
fit to the noise histogram overestimates the average and the 
number of small outliers, and it consistently underestimates the 
number of large outliers (F ig.[3ll. The skewnormal distribution 
dAzzalini & CapitaniollT999l e.g.,) 
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is more appropriate to describe the overall noise specifications. 
For the CoRoT data, we derive values of ^ = 0.14 and cj - 0.03 
for the location and scale parameters, and a value of a - 1.18, 
which determines the shape of the distribution {a - Q means 
the distribution is normal, a > means the distribution is 
right-skewed). These values imply right-skewed, leptokurtic dis- 
tributed noise, with a skewness g\ « 0.2 and an excess kurtosis 

52^0.1. 

Next, we simulate the influence of skewnormality on the pa- 
rameter and error estimation of a model 

F{ti)^ H+Aiim[2n{fti + (P)]. 

To do so, we generate three collections of light curves, each 
set containing at least 500 light curves with highly skewed noise 
{a - 10). To the first set of light curves, we add a high S/N 
monoperiodic sinusoidal signal (S/N~ 180), to the second set 
a low S/N monoperiodic sinusoidal signal (S/N~ 4), and to the 
third set a superposition of 200 sinusoids with S/N between 3 
and 190 (which mimics the CoRoT data of HD 180642). To esti- 
mate /, we use the peak frequency in the Scargle periodogram of 
each light curve. The other parameters are determined via ordi- 
nary linear regression. Subsequently, the distribution of each es- 
timator is compared with the theoretical formulae, as described 




Fig. 3. Basic properties of the noise on the data: {upper left) av- 
erage (black) and median (grey) per sample of 1000 points. A 
small temporary bump and a continuously rising trend are visi- 
ble. The fact that the median is consistently lower than the av- 
erage, suggest a nonzero skewness. {upper right) log - log plot 
of the Scargle periodogram, gray lines are straight line fits. The 
noise is mostly white, except for a low frequency steep decay due 
to correlation effects, {lower left) Histogram of the data (black) 
with a normal fit (dashed grey line) and a skewnormal fit (solid 
grey line). The skewnormal distribution fits the wings better than 
the normal distribution, {lower right) The residuals of the his- 
togram fits, after subtracting a normal fit (black) and after sub- 
tracting a skew normal fit (grey), show that a skew normal dis- 
tribution is a better estimation of the overall noise distribution 



by [Montgomery & 0'Donoghu3 ( 19991) but additionally takin g 
correlation effects into account (ISchwarzenberg-Czernvll2003h . 
e.g. the error estimate on the frequency is 
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where crr(f,) is the standard deviation of the residuals. The ef- 
fective number of observations A^eff is estimated by counting the 
average distance between sign changes in the residuals. From 
Table[ri we can see that there is no discrepancy between an es- 
timator and the real input value in the case of a monoperiodic, 
high S/N signal, besides the fact that the 'theoretical' error esti- 
mates, such as Eq. (I2]i, are rather conservative. For a multiperi- 
odic signal with low S/N, the same pattern emerges except for 
the amplitudes: the estimator of the amplitude is slightly biased 
towards smaller values, but is still well within the error bars. The 
opposite bias is found in the signals with a low S/N value, but 
here an extra bias is introduced because peaks disappear in the 
noise for low amplitude values. 

Finally, we analyze the influence of skewness on the Scargle 
periodogram in a qualitatively way using a large number of sim- 
ulations of skew normally distributed noise with different pa- 
rameters. We only find some additional noise at low frequencies, 
but this effect is only apparent for very high a values. A set of 
heteroscedastic skew normal samples also introduces additional 
noise at low frequencies, but again, the degree of heteroscedas- 
ticity has to be unrealistically high compared to the case of the 
CoRoT data, to have a significant influence. 

In conclusion, although the deviation from uncorrelated ho- 
moscedastic white Gaussian noise is significant for the CoRoT 
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data, it is not dramatic: we are dealing with a slightly right- 
skewed, leptokurtic distribution. The above simulations suggest 
that a significant influence on parameter and error estimation is 
only noticable for high departures of normality. Moreover, as 
will become clear in the following sections, the noise level is in- 
herently low compared to the analyzed signals, and we are con- 
servative in our significance criteria: in the following, we adopt 
a p value of p - 0.001 in hypothesis testing, so that an order- 
of-magnitude estimate of p is important, rather than a precise 
value. If at all, only the correlation effects are worth taking into 
ac count for our analysis; fills is done with the method outlined 
In ISchwarzenberg-Czernvl (|2003|) . 

3. Modelling of the CoRoT light curve 

Most of the calculations concerning stellar oscillations of /? Cep 
stars assume modes with small amplitudes, to be able to treat 
multiperiodicity as a linear superposition of multiple modes with 
infinite lifetime. Fitting simple sine functions, each with constant 
frequency, amplitude and phase through data represents the sim- 
plest first order deviations from a theoretical equilibrium state of 
the star However, when the perturbations are not confined to the 
linear regime, higher order effects can only be modelled when 
different sines are combined and/or harmonics are allowed for, 
spreading the signature of a nonlinear effect in a Fourier peri- 
odogram over a wide range of frequencies. 

Several physical origins of nonlinear effects in a light 
curve are plausible. These include a nonlinear response 
of the stellar flux, leading to a distortion of the light 
curve dGarrido & Rodriguezlll99 6', e.g.), nonlinear mode cou- 
pUng through resona n t interaction b etwee n different modes 
(iDziembowskil 119821: iBuchler et al.l 1 19971 e.g.), excitation 
of strang e-mode oscillations in highly nonadiabatic regimes 
jSaio et a l. 1998; Glatzel 1994, e.g.). etc. In particular, nonlin- 
ear resonant mode coupling can be distinguished from compli- 
cated beating among linear modes by checking the occurrence 
of frequency and/or phase locking, which is not expected for a 
superposition of linear modes. Nonlinear oscillation signatures 
may also include time-variable amplitudes or phases. 

Given that we are dealing with the light curve of a large- 
amplitude ySCep star, which is of unprecedented quality and 
quantity, it is not a priori clear if a linear superposition of mode 
frequencies is the best approach to treat the variability in the 
CoRoT light curve of HD 180642. Therefore, we first perform 
a traditional linear analysis of the light curve. Next, some non- 
linear models are constructed and fitted to the data, as well as 
compared with the fit assuming linear mode frequencies. This 
comparison is done by means of statistical criteria taking into 
account the number of free parameters. We thus deduce the most 
likely physical interpretation of the variability of HD 180642 
from the data point of view. 

3.1. Superposition of linear modes 

The first analysis of the CoRoT light curve of HD 180642 
was done according to the traditional method, using the linear 
Scargle periodogram (Scargle 1982) and consecutive prewhiten- 
ing, translating to a well-known model of the form 

"/ 

Fi (/,) = c + 2 A; sm[2jT(fjti + 4>j)] (3) 

j=i 

for Hf frequencies. Here, Aj, fj and (pj denote the amplitudes, 
frequencies and phases. The model was evaluated at every time 
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Fig. 4. Phase diagrams of first three independent frequencies. 
Grey lines are fits, grey circles are averages of phase bins, black 
circles are data points corresponding to successive prewhitening 
stages: (top left) five harmonics of /i - 5.4868(9)d"^ (top right) 
three harmonics of = 0.2991(7) d"', (middle left to bottom 
Wg/iO one harmonic of /3 = 6.324(6) d"', ^ = 8.409(2) d-1,/5 = 
7.254(7)6)d-',/6 = 11.811(6)d-i. 



of observation f,. At each prewhitening stage, all amplitudes, 
phases and the constant factor were refitted using the original 
light curve. This method implies a frequency resolution of the 
order of the Rayleigh limit IjT - 0.0064 d ', making the fre- 
quency determination less precise when several frequencies are 
confined to a region of this width. 

It is well known that nonlinear least squares fitting in the time 
domain while leaving the frequencies, amplitudes and phases 
free, can improve seriously the fit quality compared to the case 
where the frequency values are fixed to those resulting from the 
periodogram, but also that the success of such a procedure is 
largely dependent on the appropriate choice of good starting val- 
ues, particularly when many frequencies are present. The start- 
ing values we adopted for the amplitudes and phases are those 
that resulted from ordinary least squares regression, while we 
fixed the frequencies from the Scargle periodogram. 

It was immediately clear from the first few detected frequen- 
cies that the Scargle periodogram is not the optimal choice to 
describe or detect the variability of HD 180642, although it is 
most certainly a powerful indicator and intuitive. At least three 
harmonics of the main frequency were detected, with two more 
being marginally significant. Also the second independent fre- 
quency was best modeled with several harmonics (see upper 
right panel of Fig.|4]i. The amplitudes of the remaining frequen- 
cies are small enough to be modelled by single sines, as illus- 
trated for four of them in Fig.|4] This figure also shows that the 
light curve cannot be adequately modelled by only a few fre- 
quencies and their harmonics. 

The result of this traditional analysis is a wealth in frequen- 
cies, excluding clearly a monoperiodic model. We calculated up 
to 200 statistically significant frequencies (Table lATTT i, although 



Degroote et al: The CoRoT light curve of theySCep star HD 180642 



5 



Table 1. Comparison between errors derived using theoretical formulas with correlation correction (input), and empirically derived 
estimates of the parameters and errors (estimators), based on the results of more than 500 light curves of each 50 000 simulated data 
points. The empirically derived value for the parameters are calculated as the average outcome of the simulations, while the error is 
determined as the standard deviation. 



Set 


Type 


Frequency / 




Amplitude A 


o-(A) 


Phase (f) 




Monoperiodic High S/N 


estimator 


5.4800000 ± 0.0000004 


0.000007 


2.0000 ± 0.0002 


0.004 


+0.31703 ±0.00003 


0.0006 




input 


5.4800000 


0.000007 


2.0000 


0.004 


+0.31700 


0.002 


Monoperiodic Low S/N 


estimator 


5.47998 + 0.00002 


0.0008 


0.2153 +0.0001 


0.003 


+0.319 + 0.002 


0.06 




input 


5.48000 


0.0007 


0.0200 


0.004 


+0.317 


0.2 


Multiperiodic 


estimator 


5.4868899 ± 0.0000003 


0.000007 


34912 ± 3 


66 


-0.03552 ± 0.00003 


0.0006 




input 


5.4868900 


0.00001 


34918 


98 


-0.03551 


0.003 



it has to be noted that "only" about 100 of them would be con- 
sidered as r iot due to no ise using the traditional signal-to-noise 
criterion of iBreger et al . ( 1993). However, there is no doubt that 
the peaks are not due to noise because of two reasons. First, even 
after prewhitening 200 frequencies with a nonlinear least squares 
fit, the residual amplitudes are far above the instrumental noise 
level (discussed in Sect. l2.2l ). which would be expected if the 
signal would be composed of a superposition of linear modes. 
This can also be seen in Fig. [12] which we will discuss further 
in the text. Second, it is instructive to describe the distribution of 
frequencies across the spectrum, to see where all the frequencies 
reside. If most of the detected peaks are due to noise, they should 
be more or less randomly distributed across the analyzed fre- 
quency spectrum. In order to make the interpretation more clear, 
we decided to prewhiten a model of the first dominant mode and 
its five significant harmonics, or 

5 

F{ti) ^c + Y^aj sin[27r(j7i f,- + ^j)] , (4) 

where the initial harmonic fit was improved with a nonlinear 
least squares fit, letting also the frequency variable but with fixed 
harmonic combinations. Then, a power spectrum normalized by 
the total variance of the prewhitened data (denoted as Z) was 
calculated. This means that the expected noise level under the 
assumption of Gaussian white noise corresponds to Z = 1. The 
small deviation from this assumption (see Sect. l2.2l i. implies that 
Z = 1 slightly underestimates the true noise level. Next, the peri- 
odogram was averaged using Gaussian filters with tri = 0.1 d"' 
(to smooth out the peaks) and 1x2 = 2 d"' to estimate the empir- 
ical noise level. The result is shown in Fig.|5] Noticeable power 
excess exists around 0.3 d"', 1.0 d"', 6.3 d"', 7.3 d"', 8.4 
8.8 d-',9.8 d ', 10.4 d ', 11.0 d ', 12.3 d"\l3.9 d ', and, fi- 
nally, to a lesser extent also 14.15 d"'. Most of these power 
excess regions are not due to one large peak, but represent a 
smoothing of many closely spaced peaks in the periodogram, 
e.g., in the low frequency region (> 1 d"'). It is clear that the 
low amplitude frequencies are not due to noise, but is actual sig- 
nal; it is apparent that there are almost no frequencies nor any 
sign of power excess in the region between 1 d ' and 5 d '. The 
higher frequency regions (> 14 d ') are much closer to the the- 
oretical noise level, but are at the same time contaminated by 
secondary window peaks. 

As it turned out, several of the frequencies are linear com- 
binations of other frequencies (Table |2]i. The influence of a few 
frequencies is thus widely spread over the entire frequency spec- 
trum. This is why we consider the second model discussed be- 
low. 



3.2. Nonlinear frequency locking 

Frequencies excited through nonlinear resonant mode coupling 
can manifest themselves in a natural way through combination 
frequencies, which may seem, at first sight, independent from 
the others. Such frequency locking is one effect t hat can be de- 
rived from the Amplitude Equation formalism ([Dziembowskil 
fT982t iBuchler et alJll997t IVan Hoolst et al.lll998l e.g.l if am- 
plitudes and phases are constant in time. Following this assump- 
tion, a summary of the most obvious combination frequencies 
is given in Table|2] All combination frequencies we r e iden tified 
following the method described in iDegroote et al.l (l2009t) . We 
selected only those combinations where the difl'erence between 
the true combination value and the real value is below half of 
the Rayleigh limit Lr = 0.0064 d ^ The nonlinear leakage can 
then be viewed to spread over a wide range of the frequency 
spectrum, roughly between d"' and 20 d^' . 

This phenomenon of combination frequencie s has been de- 
tect ed previously in 13 Cep stars, e.g. in vEri (iHandler et al.l 
[2004i) and in 12 Lac (IHandler et al.ll2006l) . For these stars, only 
positive combinations were detected. It was diflicult, therefore, 
to interpret these combinations, either in terms of light curve dis- 
tortions due to nonlinear response or due to nonlinear resonant 
mode coupling. Indeed, both these phenomena would naturally 
give rise to difference combination frequencies as well as phase 
locking, which were not detected in these two stars. Moreover, 
third order combinations were not unambiguously identified, ei- 
ther because the amplitudes were too low, or they were not ex- 
cited. 

Under the assumption that the combination frequencies are 
real in HD 180642, the amplitudes, phases and independent fre- 
quencies were refitted using the previous values as starting val- 
ues, while fixing the dependent frequencies according to their 
linear combination throughout the fit: 

"/ 

Fiiti) = c +YjAk sin[27T(ftti + (l>k)] 

i=/ (5) 
+ Y^Aism[2nifiti + (pi)] 

/=! 

with 

a linear combination of two independent frequencies. It is as- 
sumed that there are n j independent frequencies, and m j depen- 
dent frequencies. 

In order to better describe the combination frequencies and 
their origin, their relative phases and amplitudes were analysed. 
Following Buchleret al. (1997) and Vuifle (2000), they are de- 
fined as 

<pr-4>c- [ni4>i + nj(f>j] 
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Frequency (d^^) 




Fig. 5. Scargle periodogram (black) after prewhitening of dominant mode and its harmonics. Overlays are Gaussian convolutions 
with (T = 0.1 d"' (light gray) and cr = 2 d"' (dark gray), which are good indicators for power excess and the empirical average noise 
level, respectively. The empirical noise level only coincides with the theoretical noise level for white Gaussian Noise at very high 
frequencies. 



and 



with the subindex / referring to the parent mode with the largest 
amplitude, j to the parent mode with the smallest amplitude and 
the index c to the daughter mode. The relative phases and ampli- 
tudes for all candidate combination frequencies, except harmon- 
ics, from Table|2]are shown in Fig.|6] 

To make the discussion more readable, we denote each 
daughter mode by a unique designation, 

dijirii, rij) 

where / and j are indices of the largest and smallest ampli- 
tude parent modes respectively (the higher this index, the lower 
the amplitude), and n,, nj are the corresponding coefficients in 
the linear combinations (Table|2|. Sum frequencies are distin- 
guished from differences by the sign of the coefficients. Several 
interesting features appear: 

1. There are six daughter frequencies which have four prop- 
erties in common: they are a sum of the dominant mode 
with another frequency, they cluster around the same rela- 
tive phase (~ 0.15), they have comparable relative ampli- 
tudes and they have the same first-order coefficients (n, - 
nj = 1): flfi,3(l, 1), flfi,4(l, 1), di,dh 1), fl?i,8(l, 1), a?i,io(l, 1) 
and iii,ii(l, 1). 

2. A similar clustering around a common relative phase is 
visible for 6 differences: ii4,8(l, -1), iii,4(-l, 1), t/i,3(2, -1), 
'^1,2(1,-1), <^5,6(2, -1) and c/i,2(2, -2), although they have 
different coefficients and different relative amplitudes. 

3. Around the second harmonic of the dominant mode, a spac- 
ing with AF = 0.29917 d"' is clearly visible. The daugh- 
ter frequencies di^i^, -1), d\;i{2, -2) have almost the same 
relative amplitudes, and are both found in two consecutive 
prewhitening stages. 

4. The daughter frequency t/i3(l - 1) has a relative phase of ;7r/2 



The result for the fit using model F2 is provided in Table |A.2l 
while a summary diagram of the independent and combination 
frequencies is shown in Fig.|2l 



Table 2. Parents {p 1 , P2) and their orders (« 1 , n2) of combination 
frequencies fc- The Rayleigh limit is equal to Lr - 0.0064 d 
Frequency values /i , /2 and fc correspond to the highest peaks in 
successive prewhitening stages. A = |ni/i +112/2- fd denotes dif- 
ference between the true linear combination and the found value. 
Column ID shows the unique designation of daughter mode /■ 
(see text). 



ID 


ni 


/i (d-') 


«2 /2 (d-') 


/. (d-') 


A(d-') 


diAh 1) 


1 


5.48689 


1 6.32482 


11.81164 


0.00006 


d,,3(-l,l) 


-1 


5.48689 


1 6.32482 


0.83794 


0.00001 


rfi,3(2, 1) 


2 


5.48689 


1 6.32482 


17.29841 


0.0002 


rf,.3(2,-l) 


2 


5.48689 


-1 6.32482 


4.64845 


0.0005 


rfl,2(l,l) 


1 


5.48689 


1 0.29917 


5.78662 


0.0006 


rfl.2(l,-l) 


1 


5.48689 


-1 0.29917 


5.18781 


0.00009 


rfi.2(2,-l) 


2 


5.48689 


-1 0.29917 


10.67458 


0.00002 


di,2(2,-2) 


2 


5.48689 


-2 0.29917 


10.37493 


0.0005 


di,2(l,-3) 


1 


5.48689 


-3 0.29917 


4.58920 


0.0002 


di,4(l,l) 


1 


5.48689 


1 8.40918 


13.89585 


0.0002 


rf,.4(-l,l) 


-1 


5.48689 


1 8.40918 


2.92159 


0.0007 


rfl,6(l,l) 


1 


5.48689 


1 6.14336 


11.63039 


0.0001 


rfl,8(l,l) 


1 


5.48689 


1 7.35867 


12.84432 


0.001 


di,io(l,l) 


1 


5.48689 


1 8.77086 


14.25740 


0.0003 




1 


5.48689 


1 6.26517 


11.75173 


0.0003 


^5,6(1,-1) 


1 


7.25476 


-1 6.14336 


1.11216 


0.0008 


^5,6(2,-1) 


2 


7.25476 


-1 6.14336 


8.36870 


0.0025 


^/4,8(1,-1) 


1 


8.40918 


-1 7.35866 


1.04985 


0.0007 


^7,0 (1,3) 


1 


7.10353 


3 0.89870 


9.79999 


0.0003 



In the case of HD 180642, we hence do see difference com- 
bination frequencies, in contr ast to the cases of th e two large- 
amplitude pCep stars vEri ([Handler et alJ l2004l) and 12 Lac 
(Handler et al. 2006). These difference frequencies are still well 
above 0.1 d ', and are thus in the regime of white noise (see 
FigOJ Such low combination frequencies are expected to occur 
with similar amplitudes as the sum combinations, for both the 
nonlinear distortion model and a nonlinear resonant mode cou- 
pling model. While we see more sum frequencies than differ- 
ences, we do reach the regime of g-mode frequencies through 
several combinations for HD 180642. We also found higher or- 
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Fig. 7. Summary of independent (full lines) and combination frequencies (dashed lines) for the nonlinear frequency locking m odel 
Fo described in the text. The frequencies detected in the ground-based photometry and spectroscopy by iBriquet et alJ (l2009h are 
indicated by closed and open triangles, respectively. 



der combinations here, up to order four (see Table|2|i, than for 
vEri and 12 Lac. 

We note from Fig.|6]that four combination frequencies have 
a much higher Ar-value than the others. This is simply due to 
the fact that these are the four combinations not involving the 
dominant mode (hence the denominator in the definition of 
is much smaller). Further, the relative amplitudes of the combi- 
nations involving the dominant mode cover a range of a factor 
ten and the phases cover the entire range [-0.5,0.5], although 
several relative phases of difference frequencies are equal within 
the error bars, and similarly for the sum frequencies. The largest 
relative amplitudes all occur for a sum frequency due to a three- 
mode resonance model involving the dominant mode. The dif- 
ference frequencies of the same three modes, if they occur, all 
have lower relative amplitude. We interpret this as due to non- 
linear resonant mode locking, as a nonlinear distortion would not 
privilege larger amplitudes for sum or difference frequencies. 

In principle, the relative amplitudes of the resonantly locked 
frequencies can help to constrain the mode degrees, because 
the geometric cancelling effect is different for different degrees. 
Unfortunately, we cannot use the relative amplitude values to 
derive the mode degrees, because all large-amplitude three mode 
resonances involve the dominant radial mode which does not im- 
ply geometric cancelling, and we have no other information on 
the degrees of the parent frequencies. 

3.3. Time dependent amplitudes and pliases 

In the previous section, we assumed the amplitudes and phases 
to be constant in time. However, the Amplitude Equation for- 
malism also allows for solutions where this is not the case 
dBuchler et al.ll997l) . Amplitude and phase modulations may oc- 
cur, which can be (multi)periodic or chaotic. The light curve of 
HD 180642 as measured with CoRoT is of such high quality, that 
it becomes possible to detect and model these variations through 
changes in the highly sampled phase profile. Although a sine 
function with five harmonics is a good fit in the phase diagram 
(Fig. Ell, it is also clear from the same figure that this model is 
only an 'average' model; in fact, the fit is not optimal to model a 
particular phase. Some phases can be modelled adequately with 
three harmonics, others need four, etc. Moreover, the minima 
and maxima seem to oscillate around an equilibrium value. 

To quantify this time-dependent behaviour, a harmonic fit 
was calculated for every covered phase of the main frequency. 



The number of harmonics to be used is determined from the 
statistic of the data with respect to the model. The number of 
significant harmonics was taken as the lowest one that achieves 
ax^ < 1-5. This number varies mainly between three and four, 
with few exceptions. 

We quantify the complexity of each phase profile by the ratio 
of the harmonic's amplitudes compared to the main amplitude. 
The higher this ratio, the more significant the specific harmonic 
is. Also, from each fit, we extract the fitted constant as an indi- 
cator for long term trends. Finally, peak-to-peak variations in the 
phases are calculated. Using these methods, we finally arrive at 
an adapted version of Eq. (HJi: 

5 

F3(ti) = c(ti) + 2 ajiti) sm[27T(jfiti + 0/f,))], (6) 

with 

c(td = C + ^ A^ sin[2;r(/,^f,- + 0^)], 

k 

ajiti) = A + ^A;'sin[27r(;/;^f,+0;')], 
<Pjitd = <D + ^A^/sin[2;r(;/^^f,+0f„')]- 

m 

For clarity, we first examine what the linear interpretation of 
this model would be. The simple model with ^ = and j, I, m - 1 
can be linearized with the assumption that A"^ < 1. Violation 
of this assumption only influences the amplitude determination. 
Linearizing (|6]l gives 

F'^iti) = A sm( Ft + <D) 

+ A" 12 sin[(F - f)ti -I- (O - 0" -H 7r/2)] 
+ A" 12 sin[(f + f)ti + (^ + <f>" - k/2)] 
+ AA't'12 sin[(f + f't')ti -H (<D -H 0'^)] 

+ AA't'12 sin[(F - f'l')ti -H (€) - / + n)] ^'^^ 
+ A^A't'lA sin[(F + f - f)ti -H (€) -H 0° - <^'^ -H 7r/2)] 
+ A^A't'lA sin[(F - /" - f)ti + - (p" - cf)''' - n/2)] 
+ A"A^IA sin[(F + f+ f't')ti -H (O -H 0" + (/>'''- n/2)]. 
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Table 3. Selection of predicted versus observed frequencies un- 
der the assumption of time dependent amplitudes and phases 
(model given by Eq. 



component 


secondary 


spurious frequency 


Error 


/i 


0.8379 


6.3248 


Observed 


/i 


-0.8379 


4.6490 


Observed 


2/i 


0.8379 


11.8168 


Observed 


2/i 


-0.8379 


10.1359 


- 


3/i 


0.8379 


17.2986 


Observed 


3/i 


-0.8379 


10.1359 


- 


component 


secondary 


spurious frequency 


Error 


fi 


1.7678 


7.2546 


Observed 


fi 


-1.7678 


3.7191 


- 


2/i 


1.7678 


12.7412 


Observed 


2/i 


-1.7678 


9.2060 




3/i 


1.7678 


18.2285 




3/i 


-1.7678 


14.6928 




component 


secondary 


spurious frequency 


Error 


/i 


2.5650 


8.0519 




/i 


-2.5650 


2.9219 


Observed 


2/i 


2.5650 


13.5388 




2/i 


-2.5650 


8.4088 


Observed 


3/i 


2.5650 


19.0257 




3/i 


-2.5650 


13.8957 


Observed 



-3.5 



4.0 



-4.5 



^ -5.0 
-5.5 
-6.0 




<pc - («,<^,-+1i|'^y)(27r/radf'^ 



0.4 



Fig. 6. {top) Relative phase as a function of frequency. Error 
bars denote 3 cr level. The parent modes are indicated by ver- 
tical lines. Horizontal lines denote a njl phase lag and lead with 
respect to the parent modes. Sum combinations are filled circles, 
difference combinations are open circles, {middle) Relative am- 
plitudes as a function of frequency, {bottom) Relative amplitude 
as a function of relative phase. 



or when f = /"^ 

F(r,) =Aisin(Ff-HOi) 

+ A2sin[(F-r)r + <l)2] 
H- A3 sin[(F -H /> H- 03] 

with, under the assumption that <s: A, 



(8) 



A^ + 



- arctan 



+ — — cos(/ 



• n/2) * A^ 



f A sin CD -H ^ sinCO + 4>" - / -H njl) 



A cos cD -H cos(<l) + I 



The results of the frequency analysis for the amplitudes, phases 
and constants for this model assumption can be found in 



Tables |A3] to lA.lOl The analysis of the first frequency of the 
time-dependent amplitude and phase for each harmonic is shown 
in Fig.[8] 

With this last expansion in mind, we can use the results of the 
frequency analysis of the amplitudes and phases to predict the 
occurrence of spurious frequencies, actually originating from the 
linear expansion of the nonlinear model, and thus not physically 
inherent to the star We performed this exercise for the primary 
component of the main radial mode and its harmonics. If we as- 
sume that the primary component f\ and its harmonics have a 
variable amplitude and phase with frequency / = 0.8379 d"\ 
then we can already explain the occurrence of several of the 
observed frequencies (see Table[3]l. Obviously, not all predicted 
frequencies are discovered. After fitting the light curve with this 
model and inspecting the periodogram, we can see that there are 
additional secondary peaks introduced which are not observed. 
When taking more frequencies into account, at least some of the 
secondary peaks seem to cancel out. This motivated us to con- 
struct several time-dependent amplitude models and to compare 
the fits. From Tables lA3] to lA. lOl we selected the most probable 
frequencies and let the amplitudes and phases vary accordingly. 
First, we tried with only one frequency. Gradually we added 
more frequencies, until we arrived at a maximum of four fre- 
quencies determining the amplitude time variability. 

The two frequencies 0.836 d ' and 1.767d ' clearly stand 
out: the amplitudes and phases oscillate on this time scale, but 
they are also found when applying model Fx given by Eq. (O, 
which means that the 'entire' light curve is also oscillating at 
the same rate. This is confirmed by analyzing the residuals after 
removing the constructed models: both of the frequencies are 
recovered. 

Despite the success of these time-dependent amplitude and 
phase models to explain several peaks in the periodogram, we 
have to compare them more rigorously with the models of the 
forms Eqs ([3]) and (|5]), which is the topic of the next section. 
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Fig. 8. Periodograms and phase diagrams for the time-variable amplitude (f^') and phase (ff') of the first four harmonics of the 
dominant mode. Because the fourth harmonic was not detected in every phase, there are less points in the lower panels for which 
the Nyquist frequency is also lower than for the other harmonics. 



3.4. Model evaluation 

To compare the different models with each other and determine 
their goodness-of-fit, we compute four evaluation statistics: the 
variance reduction (VR), both Akaike's Information Criterion 
(AIC) and the Bayesian Information Criterion (BIC) in the time 
domain, and the power reduction (PR) in the frequency domain. 
The AIC is defined as 

AIC =2^-21ni:„,x, (9) 

where X-max is a maximum likelihood estimator (MLE), n is the 
effective number of observations and k is the number of free pa- 
rameters in the model. Under the assumption of Gaussian white 
noise, we can insert the MLE of the noise variance, (f,^ = RSS/n 
with RSS the residual sum of squares. Criterion (|9]l then be- 
comes 

AIC = « ln(RSS/«) + 2k + n. 

The AIC can be calculated for every model, but is only relevant 
in comparisons; the lower the AIC, the better the model. 

Despite the fact that the AIC discourages the use of too many 
free parameters (unlike the variance and power reduction), the 
BIC is more suited when we want to stress the importance of 
simpler models over more complicated ones and thus increase 
the penalty for introducing new parameters. The BIC is defined 
as 

BIC = -21ni;niax H-'tlnn. (10) 
Analogous to the AIC, we can simplify this to 

BIC = nln (RSS/n) -Hit Inn. 

The BIC is also only relevant in comparisons, where again the 
lower the BIC, the better the model. We choose to use the BIC 



for our model selection rather than the criteria, because we want 
to favour simple physically appropriate models with the fewest 
possible degrees of freedom. The BIC is the most conservative 
criterion in this respect, because it requires that the gain in vari- 
ance reduction must be worth the cost of introducing extra pa- 
rameters. 

First, we computed the AIC and BIC for all prewhitening 
stages. Indeed, given the richness of the frequency spectrum, it 
could well be that not all variability can be modeled adequately 
with sine functions. Of course any type of variability can be con- 
sidered as such, but this may result in the use of too many param- 
eters. As can be seen on Fig.|9] the BIC sets the optimal number 
of sines to use for model fit to 127; the introduction of 3 ad- 
ditional parameters is not worth the gain in variance reduction 
any longer 

From Table|4] we can see that all statistics prefer the model 
where the combination frequencies were fixed, except for the 
variance reduction. That the latter is slightly worse is no sur- 
prise, as the number of free parameters is much smaller with 
fixed frequencies. 

All versions of the model described by Eq. (|6|l lead to a 
worse fit to the data than similar models Fi (in the sense of which 
peaks they can explain) of the form of Eq. Q. The differences 
are small in most cases, but there is a difference nonetheless. 
This was already foreseen in Table|3] because not all predicted 
peaks are detectable, hence the fit introduces spurious peaks. 
This effect can also be deduced from the PR: although more pa- 
rameters, the fits explain the periodogram less well than simpler 
models. 
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Prewtiitinmg stage 

Fig. 9. Comparison of the AIC (black) and BIC (grey) for the 
Scargle analysis of model Fi. The lower the value, the better 
the fit. The absolute values have no meaning, so the AIC and 
BIC cannot be compared with each other. Clearly the BIC gives 
a higher penalty for introducing extra parameters than the AIC. 
From frequency 127 on, the BIC discourages the use of addi- 
tional sines in the model fit (inset is a zoom). 

Table 4. Model evaluation. The AIC and BIC criteria are com- 
pared to a model with only the first frequency prewhitened. The 
power reduction with respect to the original periodogram is com- 
puted by numerical integration of the model's amplitude peri- 
odogram. Np is the number of parameters used in the model. 
The model numbers between brackets refer to the Equations in 
the text. For each alternative model Eqs ^ and Q separated 
by horizontal lines, the closest linear (original) model Eq. (O is 
computed, and fitted through nonlinear least squares. The best 
model according to the four considered statistics are emphasized 
in bold. 



tistically the best model, and also physically the more plausible 
one, we focus on the modes listed in Table|2] by performing a 
time-frequency analysis for every prewhitening stage where the 
frequency under consideration is the dominant one. 

A logical approach is to perform a wavelet analysis, adapted 
t o th e unequidistant signature of the dataset at hand (Fostgi 
1996). In a data set with a low enough noise level, one can also 
compare the shape of the detected peak p{f) in the Fourier peri- 
odogram with the theoretical shape of an infinitely stable mode 
of frequency fiy. 



Pin- 



f sin[7rr(/-/o)] 
[ 7rT(f-fo) 



with T the total time span. This method has the advantage that 
we concentrate on the most localized area possible in frequency 
space. On the downside, a peak not fitting the expected shape 
can also mean that there is a beating pattern on a longer time 
scale which is not well resolved by the data sets. A third method 
to investigate stability we applied is simply to cut the entire time 
series in half, and do an independent traditional Scargle analysis 
on both parts. 

The results for some of the frequencies in Table|2]are shown 
in Fig. (To] where the window for the wavelet transforms was 
taken as 40 days around the targeted frequencies. The dominant 
mode frequency and its first harmonic does not change during 
the entire CoRoT time series. Comparing the shapes of all the 
other peaks (some of which shown in the left panels of Fig.fTOb 
leads to the conclusion that some frequencies do not change 
their behaviour while others do. This is particularly the case 
for (i4_g(l,-l) X 1.05 d"'. The wavelet analysis hints towards 
changes in the amplitudes of many of the modes, but it does 
not allow thorough quantitative conclusions. Strong amplitude 
changes have also bee n found in the C oRoT data of the pulsat- 
ing Be star HD 49330 (iHuat et al.ll2009,) . 



Model 


AIC (%) 


BIC (%) 


VR (%) 


PR (%) 


N„ 


Eq. 15), fixed comb 


94.010 


94.065 


97.755 


25.50 


78 


Eq. ([3j, Scargle freqs 


94.023 


94.171 


97.761 


25.49 


100 


Eq. ([3j, NLLS freqs 


94.021 


94.169 


97.762 


25.49 


100 


Eq. 16), 1 freq 


95.417 


95.378 


96.873 


30.81 


40 


Eq. l3), Scargle freqs 


95.405 


95.314 


96.872 


30.87 


28 


Eq. 13), NLLS freqs 


95.404 


95.314 


96.872 


30.85 


28 


Eq. 16), 2 freqs 


95.228 


95.255 


97.019 


30.14 


58 


Eq. l3), Scargle freqs 


95.204 


95.130 


97.017 


30.24 


34 


Eq. l3), NLLS freqs 


95.204 


95.130 


97.017 


30.22 


34 


Eq. 16), 3 freqs 


94.887 


97.974 


97.255 


28.65 


76 


Eq. l3), Scargle freqs 


94.837 


94.783 


97.263 


28.84 


43 


Eq. l3), NLLS freqs 


94.837 


94.783 


97.263 


28.82 


43 


Eq. 16), 4 freqs 


94.746 


94.901 


97.355 


27.99 


94 


Eq. l3), Scargle freqs 


94.687 


94.664 


97.361 


28.29 


52 


Eq. l3), NLLS freqs 


94.687 


94.664 


97.361 


28.29 


52 



4. Time-frequency behaviour of parent and 
combination frequencies 

The enormous complexity of the power spectrum of HD 180642 
as measured with CoRoT and many of the previous remarks raise 
questions about the stability of the observed and treated modes in 
terms of amplitudes and frequencies. Nonlinear resonant mode 
coupling may give rise to variability in the frequencies and am- 
plitudes over time (iBuchler et al.lll99"7h . Given that this is sta- 



5. Long term frequency evolution of the dominant 
mode 

The Hipparcos satellite observed HD 180642 from March 12, 
1990 during almost three years. Since then, 191 observations as- 
sembled with the photomultiplier P7 attached to the 0.7m Swiss 
telescope at La Silla and to the 1.2m Mercator telescope at 
La Palma were added. Moreover, we downloaded 310 archival 
ASAS data points (Pig ulski & Poimari ski 2008). This brings the 
total time span of observations to 6814 days. The characteristics 
of the different datasets, as well as the first frequency value de- 
termined for each of these data sets are summarized in Table|5] 
and Fig.nn 

The estimator (|2| of the error on the frequency suggests that 
the precision simply scales with the total time span T, which 
would hypothetic ally lower the frequency error by several mag- 
nitudes in our case if the datasets were to be combined. However, 
we have to take into account the extremely biased distribution 
of observations in time: more than 99% of the observations are 
made during a time span only covering less than 3% of the to- 
tal time span. A natural measure for the uncertainty on the fre- 
quency is given by the width of the peak in the window function 
around f - 0. When adding only 83 datapoints to nearly half 
a million measurements, only tiny 'wobbles' are added to the 
window function and the overall shape of the main peaks is un- 
altered. 
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Fig. 10. Search for changes in the frequency or amplitude behaviour (in % of mean flux) for some of the frequencies in Table|2] Left 
panels compare the observed shape of the Scargle peak (black) with the theoretical infinitely stable peak (gray), middle panels show 
the result of a wavelet analysis, right panels show the comparison between the peak as calculated only using the first half of the 
time series (black) and the peak using only the last half (gray). From top to bottom: iii,3(l, 1) shows stability of the frequency, but 
a decrease in amplitude towards the end of the time series; fi shows signs of a decrease in both amplitude and frequency; iii,4(l, 1) 
proves to be relatively stable; <i4,8(l, -1) shows a huge rise in amplitude from virtually non existent to highly significant; t/i,2(l, 1) 
could evolve slowly in frequency. 
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Table 5. Datasets described in this paper, and frequency values for the dominant mode determined from a Scargle periodogram 
(Scargle 1982). All frequency values were also calculated via a PDM procedure (Stellingwerf 1982), a multiharmonic periodogram 
jSch warzenberg-Czernv 1996) and a nonlinear least squares fit with 5 (fixed) harmonics, and were consistent with each other within 
1 cr. By to we indicate the starting time of the observations, while T is the total timespan, n is the number of measurements, and / 
and (Tf are the frequency and error, respectively). 





to (y) 


T(d) 


n 


fid-') 


a-f{d-') 


/ i^m 


Hipparcos 


1990 


1091 


83 


5.4871 


0.00006 


63.508(2) 


Swiss 70cm/P7 


1997 


187 


20 


5.4869 


0.0002 


63.50(5) 


Mercator/P7 


2002 


1184 


171 


5.48693 


0.00001 


63.506(1) 


ASAS 


2001 


2798 


310 


5.48691 


0.00001 


63.505(9) 


CoRoT/SlSMO 


2007 


156 


379 785 


5.48689 


0.00003 


63.505(7) 
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Fig. 11. Frequency determinations of ground and space based 
observations [left to right: Hipparcos, Swiss 70cm, Mercator, 
ASAS and CoRoT). Vertical bars denote 2cr error in frequency, 
corrected for correlation. Horizontal bars denote time span for 
the frequency determination, circles denote midpoint of obser- 
vations. 



Instead of combining the datasets to model the stability of 
the frequency, we treat the datasets separately, and test if the fre- 
quency derived from Hipparcos measurements are equal to the 
one derived from the CoRoT mission, within estimated errors. 
In the CoRoT dataset correlation effects are taken into account. 
To determine how accurate the error estimation is in the case of 
the Hipparcos data, we randomly drew ~7000 samples from the 
CoRoT dataset, using the (scaled) time gaps from the Hipparcos 
measurements, to arrive at a comparable number of datapoints 
(* 100) in each sample. We derive an empirical frequency er- 
ror estimate of cry = 0.0002 d"', while formula (|2]i gives us a 
conservatively overestimated average value of cry - 0.001 d"'. 

If the frequencies for the dominant mode from Hipparcos 
and CoRoT are estimates for a common mean /, then the maxi- 
mum likelihood estimator for f is f - 5.48691 d"\ with a prob- 
ability of 0.6%. These facts suggest that the frequency of the 
dominant mode is not the same in these two data sets, but has 
instead decreased. From the datasets in Table|5] it is impossible 
to determine if the frequency change is gradual or sudden, so 
we refrain from any physical interpretation, although it would 
be naturally explained as a frequency decrease due to the ap- 
proaching of the star to the end of the core-hydrogen burning, as 
suggested by the low \ogg of 3.45 dex. 



6. Residual power 

As aheady hinted at above when discussing the three models 
Fi , Fi and Ft,, we did not yet reach the noise level when consid- 
ering all the 127 frequencies listed in Table lA. fl or the more re- 
stricted lists belonging to models F2 and Ft, given in Tables |A!9] 
to lA.lOl Figure[T2] shows the residual periodograms for each of 
the three considered models, where we took the time-variable 
amplitude and phase model F^ allowing for one frequency to 
describe this amplitude variability as this leads to the best BIC 
(TableUJi, but the result is similar for the other three cases of this 
model. As can be seen, all three models lead to residual power 
excess, although at a quite different level. Further prewhitening 
according to model F\ was done (see Table lA. ll which lists up 
to 200 frequencies) but is, according to us, not very useful for a 
physical interpretation of the f requencies. 

Recently, BeUcacem et al.l (12009.) interpreted the residual 
power of HD 1 80642 after prewhitening 9 1 frequencies with 
a model description as F\, in terms of stochastically excited 
modes due to turbulent convection. They derived a large spac- 
ing Av = 13.5;uHz or twice this value, from the residual power 
spectrum, excluding the frequencies below IBOyuHz (11.23 d"') 
and above 300yuHz (26 d ' ). We refer to their paper for a physical 
description and interpretation of such modes. 

Our model comparison shows that the residual power is quite 
different for our preferred physical model F2 than for model 
Fx considering 127 frequencies. Even though model F\ leads 
to lower residual power, it is statistically less good than model 
F2 if one takes into account the difference in degrees of free- 
dom. Also, 'natural' combination frequencies occur among the 
found frequencies, of which several are phase-locked. One may 
then wonder how frequencies which are resonantly excited and 
which may show time-dependent behaviour can be distinguished 
from stochastically excited ones, when they occur in the same 
frequency regime. In any case, the frequencies involved in phase 
locking are not expected to have a stochastic nature, as they 
would have random phase behaviour. Therefore, all the combi- 
nation frequencies whose phases are locked (6 sum and 6 differ- 
ence frequencies — see Fig.|6l) are likely not due to a stochastic 
process. They cover almost the entire range in frequency covered 
in Fig.|7] 

We estimated the large separations for the three resid- 
ual power spectra shown in Fig. [12] assuming them to be 
caused by stochastically excited oscillations (p-modes), by us- 
ing ech elle diagrams with extracted f requencies and autocorre- 
lations (IChristensen-Dalsgaardlll988h . We thus find, in the fre- 
quency range from 50 yuHz to 300 //Hz, Avi = 12. 1 ± 0.2 //Hz for 
the residuals of model Fx, Av2 = 12.9 + 0.2/zHz for the resid- 
uals of model F2 and Av3 - 18.0 + 0.2 //Hz for the residuals of 
model F3. The autocorrelation diagram for the residuals of the 
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Fig. 12. Residual periodograms for each of the three prewhitened models Fi , F2 and described in the text. 



models also shows a smaller bump around 24 //Hz and 36 //Hz. 
Assuming ^ = 1 modes as the cause of Avi , AV3 could be inter- 
preted as the distance between f = and ^ = 1 modes but we 
regard this as a tentative result which needs further confirmation. 

The detection of solar-like oscillations in this /3Cep star as 
an explanation of the residual excess power seems to be robust 
against the various models for the prewhitening of the /c-driven 
and resonantly excited modes. The value of the large spacing is, 
however, light curve model dependent. 



7. Discussion and conclusion 

The available CoRoT data of the /? Cep star HD 1 80642 provided 
us with a wealth of information. 

Beyond the previously known dominant mode, many more 
frequencies are detected. A large fraction of those does not 
change their behaviour during the time span of the CoRoT data, 
while others do. Light curve modelling using different under- 
lying functional assumptions led us to prefer a model based 
on nonlinear mode interaction, with 1 1 independent frequencies 
and 22 three-resonance combinations (among which some har- 
monics) covering the frequency range [0.3; 22] d '. This model 
selection was based on statistical criteria, without considering 
physical arguments a priori. Nine of the independent frequen- 
cies of this model are in the range expected for/? Cep stars, i.e., 
between 5 and 9 d This model is also the most logical one in 
terms of physical interpretation. Indeed, the nonlinear frequency 
locking is a natural consequence of the large amplitude of the 
dominant radial mode of the star Five of these 33 frequencies 
are in the regime of high-order g modes with frequencies below 
2 d ' for stellar models appropriate for the star The relative am- 
plitudes of the coupling frequencies differ an order of magnitude 
and seem to point towards nonlinear resonant mode excitation 
and phase locking for several of these frequencies, particularly 
for those in the g-mode frequency regime. 

Our observational results constitute a fruitful starting point 
for detailed seismic modelling of this star, particularly if some 



of the frequencies derived here could be identified. An extensive 
ground-based observing campaign has been o rganised with that 
goal and is discussed in iBriquet et al.l ( l2009h . As indicated on 
Fig.|2l it leads to fully consistent frequency results with those 
found here, with 9 high-amplitude frequencies in common. 
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Table A.l. Phase in 2:T/rad, variance reduction (VR), AIC and BIC in percentage. Amplitudes in fraction of the mean observed flux. Signal-to- 
noise (SNR) is calculated over a 6d"' in the periodogram; C = -0.00005. 
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0.000444 


2.649499 


0.005136 


0.00029 


0.00004 


-0.2926 


0. 


1259 


3, 


.7 


98, 


.70 


92.28 


93 


.25 


107 


7, 


466602 


000443 


86 419001 


005127 


00029 


00004 


-0.1099 


0. 


1257 


4, 


.0 


98, 


.71 


92.27 


93 


.24 


108 


13, 


.291190 


0.000452 


153.833214 


0.005236 


0.00028 


0.00004 


0.3146 


0. 


1283 


4, 


.9 


98, 


.71 


92.25 


93 


.24 


109 


8, 


.807314 


000449 


101 936500 


005200 


00028 


00004 


-0.4451 


0. 


1274 


4, 


.2 


98 


.72 


92.24 


93 


.24 


110 


8 


.823592 


0.000412 


102.124902 


0.004772 


0.00031 


0.00004 


-0.3130 


0. 


1170 


4, 


.5 


98 


.72 


92.22 


93 


.23 


1 1 1 


8, 


.837597 


0.000434 


102.287004 


0.005019 


0.00029 


0.00004 


-0.0539 


0. 


1230 


4, 


.3 


98 


.73 


92.21 


93, 


.23 


112 


6, 


.362414 


0.000444 


73.639052 


0.005137 


0.00028 


0.00003 


-0.1784 


0. 


1259 


3, 


.8 


98, 


.73 


92.19 


93, 


.23 


113 


8, 


.302855 


0.000444 


96.097863 


0.005135 


0.00028 


0.00003 


-0.4654 


0. 


1259 


4, 


.1 


98, 


.74 


92.18 


93, 


.23 


114 


9, 


.868206 


0.000447 


114.215342 


0.005174 


0.00027 


0.00003 


-0.4197 


0. 


1268 


4, 


.2 


98, 


.74 


92.16 


93 


.22 


115 


19, 


.382612 


0.000449 


224.335786 


0.005202 


0.00027 


0.00003 


0.0761 


0. 


1275 


7, 


.5 


98, 


.75 


92.15 


93, 


.22 


116 


8, 


.735203 


0.000451 


101.101887 


0.005221 


0.00027 


0.00003 


0.0035 


0. 


1280 


4, 


.1 


98, 


.75 


92.14 


93 


.22 


117 


0, 


.137406 


0.000446 


1.590343 


0.005163 


0.00027 


0.00003 


0.2740 


0. 


1266 


3, 


.5 


98, 


.76 


92.12 


93 


.22 


118 


8, 


.376218 


0.000444 


96.946966 


0.005136 


0.00027 


0.00003 


0.1986 


0. 


1259 


4, 


.0 


98, 


.76 


92.11 


93 


.22 


119 


1, 


.143357 


000460 


13 233304 


005327 


00027 


00003 


-0.0138 


0. 


1306 


3, 


.6 


98 


.77 


92.10 


93 


.22 


120 


10, 


.449163 


0.000467 


120.939385 


0.005406 


0.00026 


0.00003 


0.4923 


0. 


1325 


4, 


.1 


98 


.77 


92.08 


93, 


.22 


121 


7, 


.184473 


0.000455 


83.153622 


0.005271 


0.00027 


0.00003 


0.1653 


0. 


1292 


3, 


.8 


98, 


.78 


92.07 


93, 


.22 


122 


12, 


.588851 


0.000468 


145.704299 


0.005418 


0.00026 


0.00003 


0.0927 


0. 


1328 


4, 


.5 


98, 


.78 


92.06 


93 


.21 


123 


1 


.347183 


0.000473 


15.592394 


0.005479 


0.00026 


0.00003 


0.3051 


0. 


1343 


3, 


.5 


98, 


.79 


92.05 


93, 


.21 


124 


4 


375634 


000474 


50 643915 


005483 


00026 


00003 


-0.2058 


0. 


1344 


3, 


.5 


98, 


.79 


92.03 


93, 


.21 


125 


15, 


.849010 


0.000472 


183.437619 


0.005461 


0.00026 


0.00003 


0.4390 


0. 


1339 


5, 


.7 


98, 


.79 


92.02 


93 


.21 


126 


0, 


.542048 


000473 




005470 


00025 


00003 


-0.3701 


0. 


1341 


3, 


.4 


98, 


.80 


92.01 


93 


.21 


127 


0, 


.989224 


0.000468 


11.449349 


0.005411 


0.00025 


0.00003 


0.4361 


0. 


1326 


3, 


.5 


98, 


.80 


92.00 


93 


.21 


128 


4, 


.463677 


0.000470 


51.662926 


0.005437 


0.00025 


0.00003 


0.1705 


0. 


1333 


3, 


.4 


98 


.81 


91.98 


93 


.21 


129 


1 


.457477 


0.000474 


16.868952 


0.005483 


0.00025 


0.00003 


-0.2285 


0. 


1344 


3, 


.4 


98, 


.81 


91.97 


93, 


.21 


130 


0, 


.270580 


0.000468 


3.131714 


0.005412 


0.00025 


0.00003 


-0.1148 


0. 


1327 


3, 


.3 


98, 


.82 


91.96 


93, 


.21 


131 


7, 


.917871 


0.000482 


91.642025 


0.005580 


0.00025 


0.00003 


-0.2305 


0. 


1368 


3, 


.7 


98, 


.82 


91.95 


93, 


.21 


132 


1 


.296026 


0.000483 


15.000305 


0.005589 


0.00025 


0.00003 


-0.4679 


0. 


1370 


3, 


.4 


98 


.82 


91.94 


93 


.22 


133 


0, 


.237799 


0.000454 


2.752299 


0.005256 


0.00026 


0.00003 


-0.3142 


0. 


1288 


3, 


.3 


98, 


.83 


91.93 


93 


.22 
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Table A.l. continued. 



/(d-') 6^(d-') 



/ (yuHz) £f iyRz) AA/F 



SNR 



VR 



AIC 



BIC 



134 
135 
136 
137 
138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
178 
179 
180 
181 
182 
183 
184 
185 
186 
187 
188 
189 
190 
191 
192 
193 
194 
195 
196 
197 
198 
199 
200 



2, 
0, 
0. 
6, 
0. 
8. 
1, 
2, 

14. 
1, 
2, 
0. 

11. 
7. 
7. 

12. 

12. 
6, 
4, 
0, 
7. 
8. 

13. 
2. 
6, 

11. 

10. 
8, 
0, 

16. 
7. 

11. 
1. 
1. 
1, 

13. 

10. 

12. 
0. 
4. 

11. 
1. 
0. 
8, 
9, 
1, 
1, 
5, 
0, 
6. 
8. 

10. 
9. 
4, 
6, 
9, 
9, 
1, 
5. 
0. 

15. 
7. 

11. 
6, 
7, 
7. 
2. 



,043672 
,736254 
,924792 
,274253 
,453644 
,352350 
,128681 
,157032 
,257403 
,038144 
,427459 
,347300 
,751733 
,278539 
,333070 
,955309 
,844315 
,201708 
,437075 
,674955 
,295163 
,289375 
,790030 
,790366 
,330152 
,546618 
,219603 
,168456 
,658442 
161520 
,442907 
,426610 
,587122 
,767490 
,666735 
,883301 
,390805 
,741223 
,309166 
,745479 
692495 
,359556 
,975971 
,418358 
,786730 
,375485 
,177101 
,267010 
,326968 
,898981 
,457694 
,319399 
,755457 
,589204 
,132605 
,860127 
,924674 
,029558 
,751525 
,123769 
,287241 
,389585 
619288 
,229125 
,351882 
,285282 
,227274 



0.000478 
0.000447 
0.000474 
0.000455 
0.000455 
0.000463 
0.000462 
0.000474 
0.000477 
0.000470 
0.000484 
0.000475 
0.000487 
0.000482 
0.000459 
0.000487 
0.000487 
0.000483 
0.000481 
0.000478 
0.000468 
0.000476 
0.000479 
0.000486 
0.000466 
0.000482 
0.000484 
0.000483 
0.000479 
0.000492 
0.000462 
0.000493 
0.000487 
0.000488 
0.000489 
0.000486 
0.000456 
0.000485 
0.000453 
0.000489 
0.000486 
0.000488 
0.000488 
0.000461 
0.000488 
0.000476 
0.000482 
0.000497 
0.000468 
0.000499 
0.000494 
0.000499 
0.000500 
0.000504 
0.000488 
0.000485 
0.000499 
0.000447 
0.000499 
0.000475 
0.000512 
0.000473 
0.000500 
0.000508 
0.000499 
0.000476 
0.000506 



23.653608 
8.521458 
10.703616 
72.618664 
5.250507 
96.670712 
13.063439 
24.965648 

165.016241 
12.015553 
28.095586 
4.019673 

136.015424 
84.242348 
84.873501 

149.945711 

148.661057 
71.779026 
51.355038 
7.811985 
84.434760 
95.941836 

159.606828 
32.295906 
73.265645 

133.641409 

118.282442 
94.542314 
7.620860 

187.054627 
86.144753 

132.252429 
18.369468 
20.457055 
19.290919 

160.686356 

120.263944 

147.467860 
3.578306 
54.924531 

135.329801 
15.735597 
11.295957 
97.434701 

113.272343 
15.919966 
13.623856 
60.960760 
3.784349 
79.849312 
97.889973 

119.437493 

112.910384 
53.115787 
70.979230 

114.121844 

114.868909 
11.916180 
66.568575 
1.432511 

176.935655 
85.527608 

134.482497 
72.096349 
85.091228 
84.320395 
25.778639 



0.005535 
0.005169 
0.005483 
0.005271 
0.005268 
0.005356 
0.005343 
0.005490 
0.005526 
0.005439 
0.005598 
0.005495 
0.005638 
0.005576 
0.005313 
0.005633 
0.005641 
0.005588 
0.005565 
0.005533 
0.005421 
0.005509 
0.005547 
0.005625 
0.005390 
0.005580 
0.005602 
0.005595 
0.005549 
0.005696 
0.005352 
0.005705 
0.005632 
0.005645 
0.005656 
0.005626 
0.005273 
0.005615 
0.005237 
0.005663 
0.005627 
0.005650 
0.005645 
0.005331 
0.005650 
0.005515 
0.005584 
0.005752 
0.005418 
0.005781 
0.005717 
0.005773 
0.005792 
0.005829 
0.005648 
0.005615 
0.005773 
0.005176 
0.005781 
0.005500 
0.005925 
0.005471 
0.005789 
0.005883 
0.005776 
0.005512 
0.005852 



0.00025 
0.00026 
0.00024 
0.00025 
0.00025 
0.00025 
0.00025 
0.00024 
0.00024 
0.00024 
0.00023 
0.00024 
0.00023 
0.00024 
0.00025 
0.00023 
0.00023 
0.00023 
0.00023 
0.00023 
0.00024 
0.00023 
0.00023 
0.00023 
0.00023 
0.00023 
0.00022 
0.00022 
0.00023 
0.00022 
0.00023 
0.00022 
0.00022 
0.00022 
0.00022 
0.00022 
0.00023 
0.00022 
0.00023 
0.00022 
0.00022 
0.00022 
0.00022 
0.00023 
0.00021 
0.00022 
0.00021 
0.00021 
0.00022 
0.00021 
0.00021 
0.00021 
0.00021 
0.00021 
0.00021 
0.00021 
0.00020 
0.00023 
0.00020 
0.00022 
0.00020 
0.00021 
0.00020 
0.00020 
0.00020 
0.00021 
0.00020 



0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 
0.00003 



-0.0814 
0.1303 
0.3061 
-0.0798 
0.1580 
0.2023 
0.4418 
0.4815 
-0.1181 
-0.2466 
-0.1941 
-0.2060 
-0.1587 
-0.4180 
-0.1156 
0.4954 
0.1819 
-0.2326 
0.4883 
0.1507 
-0.1483 
-0.4974 
0.3534 
-0.3009 
0.0448 
0.1608 
-0.1732 
0.4513 
-0.2895 
0.4617 
-0.0362 
0.1803 
0.4358 
-0.2663 
-0.4967 
-0.0029 
0.3919 
-0.2775 
0.0626 
0.3323 
-0.0403 
0.4042 
-0.4811 
-0.2797 
-0.1297 
0.4005 
-0.2543 
-0.2691 
0.0205 
0.4291 
-0.3974 
0.3013 
-0.0400 
-0.3857 
-0.3332 
-0.0811 
0.4594 
-0.2873 
0.3417 
0.2788 
-0.1422 
0.2971 
0.4728 
0.0977 
0.2438 
-0.1116 
-0.2862 



0.1357 
0.1267 
0.1344 
0.1292 
0.1291 
0.1313 
0.1310 
0.1346 
0.1355 
0.1333 
0.1372 
0.1347 
0.1382 
0.1367 
0.1302 
0.1381 
0.1383 
0.1370 
0.1364 
0.1356 
0.1329 
0.1350 
0.1360 
0.1379 
0.1321 
0.1368 
0.1373 
0.1371 
0.1360 
0.1396 
0.1312 
0.1398 
0.1381 
0.1384 
0.1386 
0.1379 
0.1292 
0.1376 
0.1284 
0.1388 
0.1379 
0.1385 
0.1384 
0.1307 
0.1385 
0.1352 
0.1369 
0.1410 
0.1328 
0.1417 
0.1401 
0.1415 
0.1420 
0.1429 
0.1384 
0.1376 
0.1415 
0.1269 
0.1417 
0.1348 
0.1452 
0.1341 
0.1419 
0.1442 
0.1416 
0.1351 
0.1434 



3.5 
3.4 
3.4 
3.5 
3.4 
3.7 
3.4 
3.5 
4.6 
3.4 
3.4 
3.3 
4.0 
3.5 
3.5 
4.2 
4.2 
3.4 
3.3 
3.3 
3.5 
3.6 
4.3 
3.3 
3.3 
3.9 
3.7 
3.6 
3.2 
5.3 
3.4 
3.8 
3.3 
3.3 
3.3 
4.2 
3.7 
4.0 
3.2 
3.2 
3.8 
3.3 
3.2 
3.5 
3.6 
3.2 
3.2 
3.2 
3.1 
3.3 
3.4 
3.5 
3.5 
3.1 
3.2 
3.5 
3.5 
3.1 
3.1 
3.0 
4.5 
3.2 
3.6 
3.1 
3.2 
3.2 
3.1 



98.83 
98.84 
98.84 
98.84 
98.85 
98.85 
98.86 
98.86 
98.86 
98.87 
98.87 
98.87 
98.88 
98.88 
98.88 
98.89 
98.89 
98.90 
98.90 
98.90 
98.91 
98.91 
98.91 
98.92 
98.92 
98.92 
98.93 
98.93 
98.93 
98.94 
98.94 
98.94 
98.95 
98.95 
98.95 
98.96 
98.96 
98.96 
98.97 
98.97 
98.97 
98.97 
98.98 
98.98 
98.98 
98.99 
98.99 
98.99 
99.00 
99.00 
99.00 
99.00 
99.01 
99.01 
99.01 
99.02 
99.02 
99.02 
99.02 
99.03 
99.03 
99.03 
99.03 
99.04 
99.04 
99.04 
99.05 



91.92 
91.90 
91.89 
91.88 
91.87 
91.86 
91.85 
91.83 
91.82 
91.81 
91.80 
91.79 
91.78 
91.77 
91.76 
91.75 
91.74 
91.73 
91.72 
91.70 
91.69 
91.68 
91.67 
91.66 
91.65 
91.64 
91.63 
91.62 
91.61 
91.60 
91.59 
91.58 
91.57 
91.56 
91.55 
91.54 
91.53 
91.52 
91.51 
91.50 
91.49 
91.48 
91.47 
91.46 
91.45 
91.44 
91.43 
91.42 
91.41 
91.40 
91.39 
91.38 
91.38 
91.37 
91.36 
91.35 
91.34 
91.33 
91.32 
91.31 
91.31 
91.30 
91.29 
91.28 
91.27 
91.26 
91.25 



93.22 
93.22 
93.22 
93.22 
93.22 
93.22 
93.22 
93.22 
93.22 
93.22 
93.22 
93.22 
93.23 
93.23 
93.23 
93.23 
93.23 
93.23 
93.23 
93.23 
93.24 
93.24 
93.24 
93.24 
93.24 
93.24 
93.25 
93.25 
93.25 
93.25 
93.25 
93.26 
93.26 
93.26 
93.26 
93.26 
93.27 
93.27 
93.27 
93.27 
93.27 
93.28 
93.28 
93.28 
93.28 
93.29 
93.29 
93.29 
93.29 
93.30 
93.30 
93.30 
93.30 
93.31 
93.31 
93.31 
93.32 
93.32 
93.32 
93.33 
93.33 
93.34 
93.34 
93.34 
93.35 
93.35 
93.35 
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Table A.3. Variable amplitude/phase model: parameters for ay . Table A.7. Variable amplitude/phase model: parameters for (pi 
C = 0.03705 C = 0.00203 



AA/F 




f(d-') 




/ (/^Hz) 




€phi 


0.00244 


0.00011 


0.8379 


0.0002 


9.6985 


0.13 


0.05 


0.00117 


0.00011 


1.7679 


0.0003 


20.4613 


0.41 


0.10 


0.00087 


0.00011 


0.6564 


0.0004 


7.5976 


0.34 


0.10 


0.00075 


0.00011 


1.7588 


0.0005 


20.3564 


-0.39 


0.10 


0.00066 


0.00011 


1.6160 


0.0006 


18.7031 


-0.18 


0.20 


0.00064 


0.00011 


0.2999 


0.0006 


3.4708 


-0.45 


0.20 


0.00062 


0.00010 


1.8717 


0.0006 


21.6636 


-0.17 


0.20 


0.00051 


0.00010 


1.8870 


0.0007 


21.8400 


0.23 


0.20 


0.00046 


0.00010 


0.1976 


0.0007 


2.2865 


0.21 


0.20 


0.00045 


0.00010 


1.9380 


0.0007 


22.4302 


-0.41 


0.20 



Table A.4. Variable amplitude/phase model: parameters for ai- 
C = 0.00675 



AA/F 








/ (mHz) 





Ephi 


0.00131 


0.00011 


2.5651 


0.0003 


29.6883 


0.42 


0.08 


0.00085 


0.00010 


0.8377 


0.0004 


9.6958 


0.34 


0.10 


0.00075 


0.00010 


1.7678 


0.0005 


20.4611 


0.37 


0.10 


0.00061 


0.00010 


1.6166 


0.0005 


18.7106 


-0.19 


0.20 


0.00055 


0.00010 


0.5992 


0.0006 


6.9349 


0.04 


0.20 


0.00047 


0.00010 


1.8708 


0.0007 


21.6530 


-0.06 


0.20 


0.00046 


0.00008 


0.6130 


0.0007 


7.0951 


0.14 


0.20 


0.00045 


0.00008 


2.6428 


0.0007 


30.5880 


0.34 


0.20 


0.00045 


0.00008 


0.2991 


0.0007 


3.4612 


0.03 


0.20 


0.00043 


0.00008 


0.6561 


0.0007 


7.5937 


-0.44 


0.20 



Table A.5. Variable amplitude/phase model: parameters for a-}. 
C = 0.00442 



AA/F 








/ (mHz) 







0.00087 


0.00006 


2.5650 


0.0002 


29.6870 


0.46 


0.07 


0.00049 


0.00005 


0.8376 


0.0004 


9.6942 


0.22 


0.10 


0.00027 


0.00005 


0.5993 


0.0007 


6.9367 


-0.02 


0.20 


0.00026 


0.00005 


2.6425 


0.0007 


30.5840 


0.42 


0.20 


0.00026 


0.00005 


0.6569 


0.0008 


7.6031 


0.34 


0.20 


0.00026 


0.00005 


1.7681 


0.0007 


20.4638 


0.37 


0.20 


0.00024 


0.00005 


0.6126 


0.0008 


7.0899 


0.16 


0.20 


0.00024 


0.00005 


2.2028 


0.0008 


25.4959 


-0.12 


0.20 


0.00023 


0.00005 


2.3188 


0.0008 


26.8375 


-0.07 


0.20 


0.00023 


0.00005 


0.2985 


0.0008 


3.4552 


0.14 


0.20 



Table A.6. Variable amplitude/phase model: parameters for a4. 
C = 0.00203 

AA/F €A f(d-') ef /QuHz) "^T 

0.00016 0.00004 0.8375 0.0008 9.6938 -0.45 0.20 



AA/F 


ea 


/ id-') 




/ (i^m 


'P 


£phi 


0.00835 


0.00050 


0.8378 


0.0002 


9.6968 


-0.13 


0.07 


0.00539 


0.00050 


1.7679 


0.0003 


20.4616 


0.12 


0.10 


0.00482 


0.00050 


0.6564 


0.0004 


7.5976 


0.07 


0.10 


0.00436 


0.00050 


2.5644 


0.0004 


29.6802 


-0.15 


0.10 


0.00374 


0.00040 


1.6165 


0.0004 


18.7092 


-0.48 


0.10 


0.00344 


0.00040 


1.8718 


0.0004 


21.6646 


-0.44 


0.10 


0.00306 


0.00040 


1.7588 


0.0005 


20.3563 


0.30 


0.10 


0.00234 


0.00040 


1.8876 


0.0006 


21.8470 


-0.11 


0.20 


0.00200 


0.00040 


0.7411 


0.0007 


8.5776 


0.24 


0.20 


0.00197 


0.00040 


0.1982 


0.0007 


2.2940 


-0.08 


0.20 


Table A.8. Variable amplitude/phase model: parameters for ^2- 


C = 0.00203 












AA/F 


ea 


/(d-') 




f(Mm 


4> 


£phi 


0.02791 


0.00300 


0.8380 


0.0004 


9.6994 


0.11 


0.10 


0.02607 


0.00300 


2.5646 


0.0004 


29.6829 


-0.34 


0.10 


0.02034 


0.00300 


0.6572 


0.0005 


7.6063 


0.33 


0.10 


0.01615 


0.00300 


0.6111 


0.0006 


7.0729 


0.43 


0.20 


0.01317 


0.00300 


1.7676 


0.0007 


20.4588 


0.20 


0.20 


0.01133 


0.00300 


0.2992 


0.0008 


3.4624 


0.29 


0.20 


0.01123 


0.00200 


1.8724 


0.0008 


21.6711 


-0.40 


0.20 


0.01044 


0.00200 


1.1135 


0.0008 


12.8875 


0.18 


0.20 


0.01073 


0.00200 


1.7594 


0.0008 


20.3632 


0.31 


0.20 


0.01022 


0.00200 


0.5977 


0.0008 


6.9179 


0.19 


0.20 


Table A.9. Variable amplitude/phase model: parameters for (p^. 


C = 0.00203 












AA/F 


€"A 






f(Mm 





ephi 


0.03866 


0.00300 


2.5646 


0.0003 


29.6825 


-0.33 


0.08 


0.02706 


0.00300 


0.8379 


0.0004 


9.6982 


-0.01 


0.10 


0.01324 


0.00300 


0.6571 


0.0008 


7.6055 


0.25 


0.20 


0.01300 


0.00300 


1.7679 


0.0008 


20.4618 


0.28 


0.20 


0.01270 


0.00300 


1.9803 


0.0008 


22.9200 


0.17 


0.20 


0.01199 


0.00300 


0.6095 


0.0008 


7.0547 


0.50 


0.20 


0.01183 


0.00300 


2.0427 


0.0008 


23.6424 


-0.41 


0.20 


0.01106 


0.00300 


2.5777 


0.0009 


29.8348 


-0.05 


0.30 


0.01101 


0.00300 


1.1718 


0.0009 


13.5625 


0.16 


0.30 


0.01049 


0.00300 


2.1939 


0.0009 


25.3920 


-0.49 


0.30 


Table A.IO. Variable amplitude/phase model: parameters for <p4. 


C = 0.00203 












AA/F 


ea 


/(d-') 




/OuHz) 


'P 


Ephi 


0.05623 


0.00800 


0.8378 


0.0005 


9.6970 


0.16 


0.10 



Degroote et al.: The CoRoT light curve of the p Cep star HD180642, Online Material p 7 

Table A.2. Frequency table with fixed combinations between parent and daughter modes. All AICs and BICs are calculated with 
respect to the best model with only the first harmonic of the dominant mode. 



AA/F 




/O/Hz) 






eiog,„(A,) 


<l>r 




SNR 


VR 


AIC 


BIC 


0.036966 


5.486888 


63.50565 


-0.0355 


- 




- 


- 


188.7 


91.86 


100.00 


100.01 


0.006313 


10.973776 


127.01130 


-0.4967 


- 


- 


- 


- 


63.0 


94.54 


98.25 


98.20 


0.004049 


16.460665 


190.51695 


0.3490 


- 


- 


- 


- 


70.4 


95.64 


97.26 


97.18 


0.003125 


0.299178 


3.46271 


-0.3753 


- 


- 


- 


- 


27.8 


96.31 


96.54 


96.45 


0.002220 


6.324787 


73.20355 


-0.2322 


- 


- 


- 


- 


21.8 


96.64 


96.12 


96.04 


0.001568 


8.408990 


97.32628 


0.4746 


- 


- 


- 


- 


17.9 


96.81 


95.91 


95.83 


0.001441 


7.255010 


83.97002 


-0.0764 


- 


- 


- 


- 


16.5 


96.95 


95.72 


95.65 


0.001269 


21.947553 


254.02260 


0.1169 


- 


- 


- 


- 


36.6 


97.06 


95.56 


95.51 


0.001181 


11.811675 


136.70920 


-0.4118 


-4.82 


0.01 


-0.14 


0.01 


15.4 


97.15 


95.42 


95.36 


0.001003 


6.143123 


71.10096 


0.0079 


- 


- 


- 


- 


10.1 


97.22 


95.32 


95.28 


0.000999 


7.103595 


82.21753 


0.3435 


- 


- 


- 


- 


11.5 


97.29 


95.21 


95.19 


0.000924 


1.111887 


12.86906 


-0.0844 


-3.17 


0.02 


-0.06 


0.02 


9.4 


97.35 


95.11 


95.09 


0.000936 


13.895879 


160.83193 


0.2930 


-4.76 


0.02 


-0.12 


0.01 


14.9 


97.41 


95.01 


94.98 


0.000852 


7.358900 


85.17245 


0.3698 


- 


- 


- 


- 


10.2 


97.46 


94.94 


94.92 


0.000848 


1.050091 


12.15383 


0.2169 


-3.17 


0.02 


0.17 


0.02 


8.6 


97.50 


94.85 


94.84 


0.000826 


0.898758 


10.40229 


0.1134 


- 


- 


- 


- 


8.5 


97.55 


94.78 


94.77 


0.000644 


11.630011 


134.60661 


-0.1515 


-4.74 


0.02 


-0.14 


0.02 


8.9 


97.58 


94.72 


94.72 


0.000653 


0.837898 


9.69790 


-0.4485 


-5.07 


0.02 


-0.25 


0.02 


6.7 


97.61 


94.67 


94.67 


0.000579 


2.922102 


33.82063 


-0.4912 


-4.98 


0.02 


0.06 


0.02 


6.0 


97.63 


94.63 


94.62 


0.000529 


8.770800 


101.51388 


-0.0181 


- 


- 


- 


- 


6.9 


97.65 


94.60 


94.61 


0.000448 


5.187710 


60.04294 


0.4809 


-5.39 


0.02 


0.13 


0.02 


4.9 


97.66 


94.58 


94.59 


0.000426 


10.674598 


123.54859 


-0.2978 


-5.41 


0.03 


0.40 


0.02 


5.9 


97.68 


94.55 


94.56 


0.000424 


10.375420 


120.08588 


-0.2473 


-5.41 


0.03 


0.12 


0.02 


5.7 


97.69 


94.53 


94.54 


0.000345 


8.366896 


96.83908 


0.1436 


-3.60 


0.04 


0.10 


0.03 


5.5 


97.70 


94.51 


94.52 


0.000390 


6.265160 


72.51343 


0.0363 










4.7 


97.71 


94.50 


94.53 


0.000393 


17.298563 


200.21485 


0.3194 


-5.30 


0.03 


-0.36 


0.02 


9.1 


97.72 


94.48 


94.51 


0.000363 


4.648990 


53.80775 


0.1698 


-5.33 


0.03 


0.05 


0.02 


4.2 


97.73 


94.46 


94.49 


0.000339 


9.799869 


113.42441 


0.1004 


-3.36 


0.04 


0.38 


0.04 


4.9 


97.73 


94.45 


94.47 


0.000289 


5.786066 


66.96836 


0.2604 


-5.58 


0.03 


-0.37 


0.02 


3.9 


97.74 


94.44 


94.46 


0.000226 


14.257688 


165.01953 


-0.1453 


-4.91 


0.04 


-0.06 


0.03 


4.6 


97.74 


94.43 


94.45 


0.000234 


11.752049 


136.01908 


-0.1853 


-4.77 


0.05 


-0.16 


0.03 


4.0 


97.75 


94.42 


94.45 


0.000205 


12.845788 


148.67810 


0.0760 


-5.16 


0.04 


-0.17 


0.03 


4.2 


97.75 


94.42 


94.44 


0.000206 


4.589353 


53.11751 


-0.3901 


-5.72 


0.04 


-0.47 


0.03 


3.1 


97.75 


94.41 


94.44 



